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ABSTRACT 

We present deep, 170 ks, Chandra X-ray observations of Abell 2219 {z = 0.23) one of the 
hottest and most X-ray luminous clusters known, and which is experiencing a major merger 
event. We discover a ‘horseshoe’ of high temperature gas surrounding the ram-pressure- 
stripped, bright, hot. X-ray cores. We confirm an X-ray shock front located north-west of 
the X-ray centroid and along the projected merger axis. We also find a second shock front to 
the south-east of the X-ray centroid making this only the second cluster where both the shock 
and reverse shock are confirmed with X-ray temperature measurements. We also present evi¬ 
dence for a sloshing cold front in the ‘remnant tail’ of one of the sub-cluster cores. The cold 
front and north-west shock front geometrically bound the radio halo and appear to be directly 
influencing the radio properties of the cluster. 
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1 INTRODUCTION 

Merger-driven shocks, generated during the collisions of galaxy 
clusters, dissipate energy into the intra-cluster medium (ICM) 
through the injection of turbulence and acceleration of particles, 
and may therefore play an important role in the evolution of large- 
scale structure in the Universe. Merger shocks can be observed as 
temperature, pressure and entropy discontinuities in the hot X-ray 
gas and measurements of the properties of these discontinuities can 
be used to constrain transport processes in the ICM (see |Marke-| 
[vitch & Vikhlinin]|2007| for a review). However, observations of 
clear merger shock fronts, with both unambiguous temperature and 
density jumps, are rare. This is due to the requirements of a near 
plane-of-sky merger, to mitigate the effects of projection, and due 
to the inherently low X-ray surface brightness of shocked regions 
in the outskirts of clusters. To date only a handful of unambiguous 
X-ray shocks with clear temperature and density jumps have been 
discovered (e.g. the Bullet cluster, |Markevitch et al.|r2002[ Abell 
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520, |Markevitch et al.|2005[ Abell 2146, |Russell et al.|2010[ Abell 
754, [Macario et al.||20lT| Abell 2744, |Owers et al.||2011[ Abell 
2034, lowers et al.|2014| l. 

Merging galaxy clusters are often associated with large-scale 
(hundreds of kpc to a few Mpc), diffuse, steep-spectrum radio emis¬ 
sion. Two flavors of this diffuse radio emission, ‘radio halos’ and 
‘radio relics’ have been identified; the first is typically unpolarised 
to the few percent level and is distributed similarly to the cluster 
gas whilst the second is located towards the cluster periphery, of¬ 
ten elongated in morphology and more strongly polarized ~ 10-20 
per cent (see |Feretti et al.|2012| for a review). These diffuse radio 
components are synchrotron radiation from relativistic particles re¬ 
accelerated by turbulence or shocks (e.g. |Petrosian|200T| |BrunettT| 
|et al. 12004) |Brunetti & Blasi|2005| l. They are associated with the 
ICM, rather than cluster galaxies, indicating non-thermal compo¬ 
nents are present on the cluster scale. Their differing morphologies 
suggest radio halos may be associated with turbulent reaccelera¬ 
tion of particles in the cluster atmosphere, while the location of 
radio relics and their relation to surface brightness edges in X-ray 
observations suggest they trace shock fronts generated by major or 
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Figure 1. Abell 2219: A three colour Subaru optical image using the Rc, V and B filters and the Chandra X-ray surface brightness image (purple). 


minor mergers (see review by |Bruggen, van Weeren & Rottgering] 
|2011| and references therein). 

In this paper we study one such merging cluster in detail, Abell 
2219 {z = 0.2256), one of the hottest (Tx — 9.5 keV, 


Allen 


& Fabian|1998 ► and most X-ray luminous clusters known ( Ebel 


ing et al.||1998L Abell 2219 is more than twice as bright in the 


X-rays as either the Bullet cluster or Abell 520. X-ray and optical 
observations show clear evidence for a disturbed, elongated mor¬ 
phology and more than one dominant mass clump (|Allen | 


Bezecourt et al.||2000l [Boschin et al.||2004l [Smith et^ 


1998 


2005 


der Linden et al.|2014| the optical and X-ray morphology of 


the Abell 2219 system are shown in Fig. [^. From a multiwave¬ 
length analysis, [Boschin et al.j ( |2004| ) suggest a complex merger 
history in Abell 2219, with the system undergoing an in-fall of 
many clumps aligned with a filament in the foreground, oriented 
at an oblique angle with respect to the line-of-sight. Abell 2219 
also hosts a well known radio halo ( [Giovannini, Tordi & Feretti] 

|1999[[Bacchi et al.|2003[|0rru et al.|2007| ) and three strong radio 

galaxies |Owen, White & Burns|1992|t. Using the Chandra X-ray 
space telescope. Million & Allenj l 2009 j l detected a large scale and 
very high-temperature shock front ~2 arcmin (~430 kpc) from the 
cluster core. Their observed temperature and density jumps lead 
to Mach numbers of >1.3 with an estimated shock velocity of 
~2500 kms“^. Surprisingly, the authors found Abell 2219 also 
appears to contain a ‘hot core’ (>20 keV). The astrophysical cause 
for this hot core is unclear. 


Here we present an analysis of new deep ~170 ks Chandra 
observations of Abell 2219. The Chandra observations and data 


reduction are discussed in Section Results from imaging, ther¬ 
modynamic maps and spectral and surface brightness profiles are 
presented in Sections [^|^and[^respectively. We discuss the X-ray 
results in the context of other multi-wavelength data on Abell 2219 
in Section [53] andand summarize our findings in Sectionj^ 

We assume a standard ACDM cosmology with Ho = 71 km s^, 
Ctrn = 0.27 and Qa = 0.73. For this cosmology, at the redshift of 
Abell 2219, an angular size of 1” corresponds to a distance of 3.588 
kpc. 


2 CHANDRA DATA PREPARATION 

Deep Chandra observations of Abell 2219 were made between 
2012 May 26th and 2012 October 15th. The observations were 
made with the Advanced CCD Imaging Spectrometer (ACIS) in 
VFAINT mode and using CCDs 1-4 and 6. The data were repro¬ 
cessed using Cl AO 4.5 and CALDB 4.5.9. Existing ACIS-S data 
taken on 2000 March 31st was also reprocessed. The total cleaned 
exposure time of the observations is 169 ks. A description of the 
observations is given in Table 

We extract and filter the lightcurves from all source-free 
CCDs, using energy ranges of 2.5 to 7 keV and 1024s bins for the 
back illuminated chips and 0.3 to 12 keV in 256s bins for front il¬ 
luminated chips, respectively, also excluding point sources. These 
binnings match the cleaning set-up used for the blank sky back¬ 
ground files. We use a sigma clipping algorithm to search for and 
reject fiares. This is followed by a visual inspection of the light 
curves. 


© 0000 RAS, MNRAS 000, 000-000 

































Shocks and edges in A2219 3 


Table 1. Description of observations. OBS ID 7892 has a cleaned exposure 
time of less than 5 ks and as such was not included in the analysis presented 
in this paper. 


OBS ID 

Date 

Instrumental 

setup 

Cleaned 
exposure (s) 

896 

2000-03-31 

ACIS-S, 235678 

40,250 

13988 

2012-05-26 

ACIS-I, 01236 

8,670 

14355 

2012-06-28 

ACIS-I, 01236 

26,100 

14356 

2012-10-15 

ACIS-I, 01236 

41,450 

14431 

2012-05-27 

ACIS-I, 01236 

34,700 

14451 

2012-06-26 

ACIS-I, 01236 

17,660 


Blank-sky backgrounds are processed and re-projected to 
match the observations. The blank-sky exposures were normalised 
to match the particle background-dominated count rate in our cor¬ 
responding target exposures, using the data in the 9.5 to 12 keV 
band. 

Narrow- and broad-band images between 0.6 and 7 keV were 
created, background subtracted and exposure corrected. Exposure 
map weighted PSF images were examined with the CIAO routine 
WAVDETECT to identify point sources. A manual inspection of the 
point sources was then performed and these sources were excluded 
from the remainder of our analysis. 

We check the blank sky backgrounds against source free re¬ 
gions of our observations for every OBS ID. We find that on the 
S5 chip of the 896 ACIS-S data an additional soft component is 
required to model the background. Fitting an APEC ( [Smith et~^ 
[200 Ij ) model for thermal Bremsstrahlung and line emission to this 
soft component, we find a best fit temperature of 1.1 keV. Over¬ 
laying ROSAT X-ray data and Digital Sky Survey optical images 
shows that there is an excess of both X-ray emission and a pro¬ 
jected galaxy over-density North (N) of the S5 chip, indicating the 
soft excess may be emission from a group. When the region of the 
chip containing most of the ‘group’ emission is excluded, no soft 
excess is required to fit the data. 

We also examined the effect of adjusting the applied back¬ 
ground normalisation by ±5 per cent and inspect the effect on the 
resulting thermodynamic maps. Within their uncertainties our mea¬ 
surements are not affected. 


3 IMAGING ANALYSIS 

The X-ray emission is elongated in the north-west (NW) to south¬ 
east (SE) direction (Fig. [^. Two bright concentrations of galax¬ 
ies are observed along the same axis, which is the likely projected 
merger axis. 

The mosaiced, exposure-corrected X-ray image, with a 2D el¬ 
liptical beta model fitted and subtracted (double beta models are 
not formally better fits), is shown in the top panel of Fig. This 
reveals two surface brightness edges (marked in cyan on the image) 
along the projected merger axis on either side of the centroid. The 
centroid of the 2D beta model is SE of the brightest region of X-ray 
emission between the two dominant galaxies and is marked with an 
‘X’. 

The bottom panel of Fig. shows an unsharp-masked X-ray 
image (made by subtracting an image smoothed by a Gaussian with 
cr = 35 from an image smoothed by a Gaussian with cr = 9) of the 
core, overlaid on a three colour optical Subaru Suprime-Cam im¬ 
age. A complicated structure is revealed in the core with two bright 



Figure 2. Top: Mosaiced, exposure corrected X-ray image, with a 2D el¬ 
liptical beta model fit and subtracted. Two surface brightness edges are ob¬ 
served and are indicated with cyan arcs. The dashed box is 3 arcmin by 
3 arcmin and is the same size as the lower panel. Bottom: Colour com¬ 
posite of an unsharp-masked X-ray image (blue) (made by subtracting an 
image smoothed by a Gaussian with cr = 35 from an image smoothed by a 
Gaussian with a = 9) and a three colour optical Subaru image. Two X-ray 
bright regions are observed in the core and are associated with the dominant 
galaxies, with the ‘peak’ in the X-ray brightness associated with BCGl. The 
positions of the dominant galaxies and the centroid of the 2D beta model 
are indicated. 


regions of X-ray emission associated with both the brightest clus¬ 
ter galaxy (BCG 1) and the second largest galaxy (BCG 2) of the 
A2219 system. These galaxies have a projected separation of ap¬ 
proximately 50 arcseconds (~180 kpc). |Boschin et al.] ( |2004| ) find 
a line-of-sight velocity difference between BCG 1 and 2 of 1,320 
kms“^. 
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Figure 3. From top left, clockwise: The temperature (keV), pseudo density (cm“^ pseudo pressure (keV cm“^ and pseudo entropy (keV 

cm2 /i/3^ maps of Abell 2219, out to 200 arcseconds (~720 kpc). These figures display 3,600 counts per bin, which at these temperatures allows us to place 
<15 per cent uncertainties in the temperature. Relative abundances are fixed to 0.3 times the solar value. The maps are the same scale, with N up and E to the 
left. The positions of the BCGs are marked with blue crosses on the temperature map and the positions of the surface brightness edges discussed in Section 
l^are marked with cyan arcs. A further density discontinuity (see Fig.|^ is marked in red. Montecarlo maps showing the uncertainty and covarience in the 
temperature and density maps are located here: http://www.web.stanford.edu/~rcanning/a2219.html 


4 THERMODYNAMIC MAPS 

In Fig. we show projected maps of temperature (in units of keV), 
pseudo density (cm“^ pseudo pressure (keV cm“^ 

and pseudo entropy (keV cm^ for Abell 2219 (where I is nor¬ 
malised to 1 Mpc). The CONTBIN routine ( |Sanders|2006| ) was used 
to bin the 0.6-7.0 keV images, with point sources masked, to a 
signal-to-noise of 60, requiring > 3,600 background-subtracted 
counts per bin. The spectra were fitted between 0.6-7.0 ke\Qwith 
an APEC thermal model ( [Smith et al.|2001|) , with a Galactic equiv- 
alent column density of Nh — 1.76 x10^^ cm“^ |Kalberla et al.j 


^ Higher energy cut-offs were also used but did not change the results. 


|2005j ). At the high systemic temperature of Abell 2219 (~ 10 keV) 
this binning allows us to place <15 per cent uncertainties on the 
projected temperature measurements and < 5 per cent uncertainties 
on the pseudo density measurements in each bin. 

The metallicity was fitted using larger regions with 20,000 
counts per bin. However, due to the high cluster temperatures, 
20,000 counts per bin only allows for 30 per cent statistical errors 
on the abundance. Within the uncertainties, the fitted abundance 
values are consistent with a constant value of 0.3 across the sys¬ 
tem, with only a slight 2a) indication for enhanced metallicity 
at the position of BCGl (and the X-ray ‘peak’). The maps shown 
in Fig.j^have relative abundance fixed to 0.3 times the solar value. 

The top left hand plot of Fig. [^reveals a ‘horseshoe’ of high 
temperature gas, which surrounds the core, and is oriented towards 
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Figure 4. The xspec normalisation and temperature in 5 equal sized re¬ 
gions spanning the core. The regions are shown in the top panel. The second 
and fourth regions coincide with the second and first BCGs, respectively. 
The density increases rapidly to a peak on the larger, dominant BCG from 
the NW, then declines slowly to the SE due to the second ‘peak’ from the 
smaller cluster. The gas temperature is high at the dominant BCG (BCGl), 
while at the second BCG the temperature is consistent with the neighbour¬ 
ing regions within the statistical uncertainties. 

the NW. This feature, with a temperature >14.5 keV, extends over 
8 independent regions in the map. In the NW, the projected location 
of the high temperature gas is roughly coincident with the NW sur¬ 
face brightness edge discussed in Section[^(marked in cyan on Fig. 
^ and coincides with a the shock front first reported by |Million &| 
Allen| ( [200^ . The entropy and density maps both have relatively 
sharp features at this location, with the entropy rising just beyond 
the edge. SE of the core, a tail of low temperature gas is apparent. 
The tail is coincident with low entropy, low pressure gas. 

The entropy map is asymmetric: it has a highly elongated 
structure, orientated NW-SE, that rises more quickly in the NW. 
The pressure map is ‘rounder’ although the pressure gradient drops 
more rapidly towards the SE. Regions of high pressure are observed 
perpendicular to the projected merger axis. 

The core is complex, with a higher density and lower entropy 
than the surrounding gas but also with a high temperature and pres¬ 
sure. The positions of the two dominant galaxies, are indicated on 
the temperature map with blue crosses. They both coincide with 
high temperature regions. To establish whether these high tempera¬ 
ture peaks are significant, we extract point source removed spectra 
from 5 regions, of equal size, spanning the core and orientated par¬ 
allel to the projected merger axis (see Eig. |^. We cannot provide 
strong constraints on the metallicities or any multiple temperature 
structure of the gas. However, single temperature APEC fits show 
that the surface brightness peak, coincident with BCG 1, is signifi¬ 
cantly hotter and denser than its surroundings. The density declines 
more slowly towards the SE of BCG 1 due to the contribution as¬ 
sociated with BCG 2 and the second subcluster. At BCG 2, the 
temperature is high but within the statistical uncertainties of the 
neighbouring regions. 


5 PROFILES 

The imaging analysis and thermodynamic maps show two clear sur¬ 
face brightness edges at large radii (> 400 kpc) along the direction 
of the projected merger axis. There are also complex, high pressure. 


high temperature regions surrounding the core. We examine these 
inner and outer features in more detail in this section. 

Whilst the intermediate, 1.8-3.0 keV, band provides good den¬ 
sity sensitivity, independent of temperature, for high temperature 
gas, the 0.6-7.0 keV full-band image is also relatively insensitive 
to variations in the gas temperature. Due to the improvements in 
signal gained by using the full band we therefore choose to make 
our surface brightness profiles using the full band. The counts are 
extracted separately on the merged counts image and merged back¬ 
ground image. These are divided by their respective exposure maps 
and the background counts are subtracted. Uncertainties are propa¬ 
gated in quadrature. 

Aside from where specifically stated, the thermodynamic pro¬ 
files presented are all deprojected profiles that assume spherical 
symmetry. We do not combine the spectra across different chips; 
instead we fit the data from each OBS ID and chip simultaneously 
in XSPEC. 


5.1 The core 

Eig. shows the surface brightness profiles of 8 regions surround¬ 
ing the core, extracted in 2 arcsecond radial bins. The centre and 
position angle of the elliptical profiles are those of the fitted cen¬ 
troid of the best fit elliptical 2D beta model discussed in Section 

m 

The steepest surface brightness profiles are in the N and NE, 
offset slightly from the direction of the projected merger axis. SE 
of the core, the surface brightness exhibits a distinctly shallower 
profile. A series of edges surround the core and correspond well 
with regions of high projected pseudo-pressure and high temper¬ 
ature. These edges are shown on the profiles in Fig. indicated 
by dashed red lines, and their positions in 2D are indicated on the 
projected temperature map, in the center of the figure. 

The surface brightness edges in the core are sharp, narrow fea¬ 
tures. It is therefore difficult to examine the spectral properties of 
these regions. The requirement of ~5000 counts to measure the 
deprojected temperature limit us to using relatively large regions 
from which to extract spectra. Fig. shows the spectral properties 
extracted in regions towards the W of the core where the highest 
projected ‘inner’ temperatures and pressures are observed. At least 
three consecutive surface brightness edges are observed in the same 
direction at radii of ~25, ~40, and ~55 arcseconds (see Fig. 1^. An 
increase in both projected and deprojected temperature is observed 
at ~30 arcseconds (~108 kpc) from the core. Discontinuities or 
edges in the surface brightness, density, pressure and entropy are 
also observed, indicative of a shock. 


5.2 The outer edges 

Initially, we select by eye sectors which trace the edges observed 
in Fig.|^from which to extract surface brightness profiles and map 
the spectral properties. The surface brightness and thermodynamic 
profiles in these sectors are shown in Fig. in regions of 5,000, 
8,000 and 10,000 counts per bin. 

The thermodynamic profiles in the NW show a significant 
temperature jump 100-140 arcseconds from the core. The tempera¬ 
ture discontinuity is accompanied by a density, entropy and pres¬ 
sure discontinuity. The direction of the temperature and density 
jumps indicate a shock front. 

In the SE, a possible temperature increase is observed at ~50 
arcseconds in close proximity to a surface brightness edge and the 
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Figure 5. Surface brightness profiles extending to 100 arcseconds (~350 kpc) elliptical radius from the X-ray centroid of Abell 2219. The figure is orientated 
such that projected merger axis (NW to SE) is vertical with NW being up. The core surface brightness profiles are complex with many edges observed. The 
projected temperature map is shown in the centre and many surface brightness edges correlate with regions of projected high temperature. To guide the eye, we 
have marked the obvious surface brightness edges with red dashed lines on the profiles and the corresponding arc on the temperature map. The centre, major 
and minor axis ratio and position angle are from the best fit elliptical beta model. This is centred at RA 16^40’20.173’ and Dec +46^42’30.60” with position 
angle 6> = 40.8° W from N. 
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Figure 6. Top; Temperature map of the core of Abell 2219. There are a 
minimum of 3,600 counts per bin. Bottom: Deprojected temperature, den¬ 
sity, pressure and entropy profiles in the core of Abell 2219. The red points 
are extracted from the annular regions shown in the top panel (red), which 
have a minimum of 8,000 counts per region. The black points are extracted 
from the same sectors but with 5,000 counts per region. The green points 
show the projected temperatures in regions with 8000 counts. All abun¬ 
dances are fixed at 0.3 solar. There is weak evidence for temperature and 
density discontinuities in an arc roughly ~20-30 arcseconds from the core, 
which coincides with the regions of high projected temperature in the map 
above. Three surface brightness edges exist at ~20, 40 and 55 arcseconds. 


location of BCG 2 (dashed line). A third possible density disconti¬ 
nuity is also apparent at a radius of ~ 150 arcseconds. 

We further refine our surface brightness extraction regions by 
dividing the NW and SE sectors into three, measuring the strength 
of the surface brightness ‘edge’ in each of the three sectors and 
varying the position angles over which they are extracted. Though 
we have tried to mitigate them, uncertainties in estimating the el- 
lipticity of the edge and the centring of the extraction regions will 
exist. These effects blue the edges leading to less contrast in surface 
brightness and thermodynamic features (unfortunately our signal- 
to-noise is not high enough at the edge to bring an edge model 
directly to our X-ray image and fit for these parameters). Our 
‘sharpest’ surface brightness profiles are shown in Fig.[^(the width 
of the annuli are chosen to be a minimum of 4 arcsec or alterna¬ 
tively a minimum of 100 counts per region). 

We extract spectra from the same sectors with the requirement 
that each of the annuli from which spectra were extracted must have 
a minimum of 2,000 counts per region in the 0.6-7.0 keV band 
(calculated from the background subtracted, point source removed 
images). Since, the temperature is not well constrained with only 
2,000 counts, we tie the temperature in adjacent regions to a mini¬ 
mum of 6,000 (black points) counts per region. We use the surface 
brightness profiles to ensure the temperatures are not tied across a 
surface brightness edg^ The profiles are shown in Fig.[^ 

The shock front is clearly seen towards the NW (~ 150”) and, 
surprisingly, two density discontinuities are seen in the SE. The 
farthest (~ 170”) corresponds to the edge observed in the unsharp- 
masked image (Fig.[^, but an additional discontinuity, not seen in 
the surface brightness profile, is also observed at ~ 100 arcseconds. 

To estimate the edge properties, we model the density on ei¬ 
ther side of the discontinuities as a power-law of the form: 

( 1 ) 

where rii is the density, ai is the power law index, Ni is the normal¬ 
isation at each side, z, of the discontinuity and 77/ is the position 
of the discontinuity. This can be fit directly to the deprojected den¬ 
sity profile or projected and fit to the surface brightness profile (e.g. 
[Owers et al.|2009| ). To model the surface brightness we make the 
assumptions of spherical symmetry, that the underlying beta model 
can be neglected and that the surface brightness is not sensitive to 
the temperature variations about the shock (as mentioned in Sec¬ 
tion]^ at the high temperatures of Abell 2219 and in the X-ray 
energy-band considered, this is a reasonable assumption). 

The best fits to the surface brightness profiles and to the de¬ 
projected density profiles are shown by the red lines in Fig.[^ The 
uncertainties on the fits are found using an monte-carlo approach 
and the errors quoted on Fig. are 95%. However, it should be 
noted that changing the region over which the fit is performed in¬ 
fluences the best fit. In practice, we must use regions which are 
large enough to have many data points to fit to either side of the 
jump but small enough that the assumption of power-laws in den¬ 
sity is valid and that other sharp features in the surface brightness 
profiles are not included. The regions used for the fits can be seen 
on Fig. To estimate the systematic uncertainty we perform the 

^ Deprojecting the profiles in low surface brightness regions can introduce 
oscillatory noise into the neighbouring bins. To mitigate this we used the 
hard band (4.0-7.0 keV) image to determine the radii to which we should 
extend our measurements. A 12 keV APEC model at the redshift of our 
cluster has 1.4 per cent of the counts in this energy range compared with 
the 0.6-7.0 keV band. We therefore require our profiles to have >50 counts 
in the hard band 
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Figure 7. The surface brightness, deprojected temperature, density, pressure and entropy profiles out to 200 arcseconds. Regions with >8000 counts are 
indicated in blue and red on the top image and in the bottom profiles. The black points have >5,000 counts per region, and the magenta and cyan diamonds 
have >10,000 counts per region. The green central temperature and the green points in the temperature profiles are projected temperatures from the 8,000 
count annular regions. The abundances are fixed at 0.3 solar. A series of edges discussed in Fig.[^can be seen within 50 arcseconds in the surface brightness 
profiles and a clear temperature jump indicating a shock front, likely the bow shock, is observed at >100 arcseconds in the NW. The SE is cooler and the 
profiles are smoother. The position of BCG 1 and BCG 2 are indicated with dashed lines in the lower panels. 
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Figure 8. Surface brightness and thermodynamic profile through sectors optimised to enharnce the surface brightness edges (error bars are 90%). The left 
hand panel shows the profiles across the SE edge and the right hand panel across the NW edge. Top: The X-ray surface brightness profiles and a residual 
profile after subtraction of a single ^-model fit. The residuals are shown here as deviations from the smooth /5-model, indicating the location of the edges is 
more easily observed in the residual profile. Middle: Deprojected density profile (2,000 counts per region). Two discontinuities are observed in the SE. While 
not clear from the surface brightness profile itself the discontinuity does correspond to a deviation from the smooth profile observed in the residual profile. 
Bottom: Deprojected temperature profile (black - 6,000 counts per region; red - 12,000 counts per region). Fits to the density discontinuities and surface 
brightness edges are shown with solid red lines. 


fit over wider and narrower regions (3 data points either way); the However, we cannot rule out a smaller temperature jump consistent 

systematic uncertainties on the Mach numbers are typically 0.15. with M = 1.2. 

The results from these fits indicate: 


(i) SE 1 shock front (~155”, 556 kpc): A fit to the surface 
brightness edge gives M = 1.16 whilst the density discontinu¬ 
ity gives M = 1.37. Taking into account the systematic uncer¬ 
tainties the Mach numbers are consistent. A simple F-test rules 
out a straight line fit to the density discontinuity with 99% prob¬ 
ability, and AT = ^ at the interface is also inconsistent with 
unity of temperature to greater than 0.01%. The best fitting value 
of AT = 0.7, indicates a Mach number of M = 1.2, and the sign 
indicates a shock front. A Mach number of 1.37 from the density 
discontinuity leads to a velocity of the shock of ~2,350 kms“^ 
for a T ~ 10 keV gas. The lower Mach number of 1.16 reduces the 
velocity to ~2,000 kmps. 

(ii) SE 2 cold front (~100”, 359 kpc): The projected density 
model gives a reasonable fit to the surface brightness, but, a simple 
F-test cannot reject the null hypothesis that no discontinuity in den¬ 
sity exists or that no edge exists in the surface brightness profile. 
However, the temperature measurements across the discontinuity 
are inconsistent at the 3% level, with the best fit AT = 1.15. 

(iii) NW shock front: (~166”, 595 kpc) The density discon¬ 
tinuity is inconsistent with no jump at the 99% level. A fit to the 
density discontinuity gives M = 1.22, consistent with the sur¬ 
face brightness edge. This leads to a velocity of the shock for a 
T ~ 10 keV gas of ~2100 kms“^. The projected distance from 
the centroid is 166 arcsec (595 kpc), which allows us to estimate 
the time since core passage of the shock as 0.26 Gyrs. The best fit 
value of AT = 0.63, suggests a higher Mach number of 1.4, the 
temperatures are inconsistent with no variation at the 0.01% level. 


5.3 Radio-X-ray connection 

Fig. shows the morphology of the diffuse (top panel, 325 MHz 
VLA radio contours reproduced from |Orru et al.|2007| l and point 
source (middle panel, 1.4 GHz VLA radio contours) radio emis¬ 
sion overlaid on the X-ray projected temperature map. The loca¬ 
tion of the fronts are labelled. Interestingly, the radio halo appears 
to be bounded by the shock front in the NW and the cold front in 
the SE (Fig.[^. Merger radio halos are thought to be generated by 
turbulence (e.g. |Brunetti et al.|2001[ Petrosian|200T]) while merger 
radio relics, or radio ‘gischt’ ( |Kempner et al.||2004| t, are located 
at the periphery of clusters and associated with particle accelera¬ 
tion at shock fronts (e.g. [Ensslin et al.||1998| l. Despite these defi¬ 
nitions, there are examples of radio halo edges coincident with X- 
ray shocks (e.g. Coma cluster, [Brown & Rudnick|20TT Abell 520, 


Markevitch et al.|20^ Abell 665, |Govoni et al.|2004[ |P eretti et ah] 


2004 [ Bullet Cluster, [Markevitch et al. 12002)[Shim well et al.|2014| l. 


5. 3.1 NW and SE X-ray shock fronts 

The morphology of a radio surface brightness edge corresponds 
well to the high temperature X-ray gas at the NW shock front. How¬ 
ever, while the high temperature gas has a clear horseshoe shape the 
radio emission has a more limited extent. These differing extents 
may be a reflection of the increased gas density along the projected 
merger axes (Eig. [^, the distribution of the pre-existing electron 
populations or an indication of where turbulence is most efficiently 
generated behind the passing shock. 
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Figure 9. X-ray projected temperature map overlaid with contours' of dif¬ 
fuse radio emission (325 MHz radio brightness contours reproduced from 
|Orru et al.|2007) (top) and VLA 1.4 GHz emission indicating the radio 
galaxies (middle). Bottom: The HST optical imaging on the same scale. 
The position of the ‘outer’ edges, discussed in Section [5^ are indicated by 
a red arc (cold front) and cyan arcs (shock fronts). 


The surface brightness and density discontinuities of the NW 
shock are consistent with a low Mach number of ~1.2 while the 
temperature favours M ~ 1.4. If the edges suffer from projection 
effects the density Mach numbers could be underestimated (e.g. 
[Markevitch, Vikhlinin & Mazzotta|2001[|Ogrean et al.|2013| >. The 

SE reverse shock also exhibits a low Mach number of 1.1 — 1.4, 
with the temperature jump in this case providing a lower Mach 
number than the density discontinuity. However, the temperature 
measurement is likely contaminated by the colder tail of stripped 
gas emanating from the disrupted sub-cluster core. Unlike the NW 
shock front the SE shock front is not associated with strong diffuse 
radio emission. 

The morphological connection between the radio halo and 
NW shock front in A2219 is strengthened by the observation that 
the spectral index steepens, from ~0.8, in the halo, to ~ 1.4, at the 
location of the front (see |Orru et al.|2007] their Eig. 9). If we make 
the assumption that the radio emission near the shock is acceler¬ 
ated by Eermi acceleration of the shock front the expected radio 
spectral index is a = ((2(M^ + 1)/(M^ — 1)) — l)/2 + 1/2, 
where M is the shock Mach number ( [Blandford & Eichler|1987| ). 
Eor a spectral index of 1.4, we require Mach numbers >2.5, sig¬ 
nificantly higher than measured. The density of the shocked gas 
would need to be 2.7 times the density in the post shock region 
whereas we measure 1.5 (a density ratio of this amount could oc¬ 
cur if the post shock gas takes less than 20% of the volume of our 
measured bin). This argues against direct acceleration, as does ar¬ 
guments about the radiative lifetime of the electrons ( 10® yrs) 

which is several times shorter than the shock crossing time in this 
system (e.g. |Brunetti|2003| l. Moreover, we see that the radio emis¬ 
sion extends beyond the shock front. Re-acceleration of a existing 
population of energetic electrons may be occurring: these electrons 
could originate from an existing radio halo or have been generated 
by the three known radio galaxies ( [Owen, White & Burns|1992] l. 
The lack of radio emission associated with the SE shock may indi¬ 
cate that turbulence is low behind the in-falling subcluster and not 
yet widespread throughout the merging system. 

The spectral index map of the radio emission appears not to 
be smooth, with regions N and W of the core steeper than the S 
(localised variations in the spectral index of the Bullet Cluster were 
also observed by |Shimwell et al.|2014| l. This may indicate that the 
halo is still forming and the diffuse radio emission is being varied 
locally by the proximity to shocks and the amount of turbulence. N 
and W of the core of A2219 we detect multiple edges with the char¬ 
acteristics of shocks (see Section [6T| . Although these have very 
low measured Mach numbers they may influencing the radio emis¬ 
sion in these central regions. 

5.3.2 SE X-ray cold front 

In the SE, the radio edge coincides very well with the SE cold front 
providing further evidence that this is a real feature. The radio emis¬ 
sion sits inside the edge, as is observed in ‘sloshing’ cold fronts in 
relatively relaxed systems (e.g. |Mazzotta & Giacintucci|2008| l. The 
abrupt edge to the halo supports the hypothesis that this ‘remnant 
tail’ gas is not yet mixed with the ICM. The spectral index map of 
|Orru et ak] ( |2007| ) suggests that the spectral index steepens towards 
the SE boundary (~ 1 zb 0.15) where we find the X-ray cold front. 
We note here the uncertainties on this steepening are large. 

In merger cold fronts the radio emission is amplified outside of 
the front, where the magnetic field is being swept up and stretched 
out, by the motion of the front through the hotter and less dense 
medium ( |Lyutikov|2006| ). In cold fronts generated by sloshing of 
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dense X-ray cores the shearing motions both amplify the mag¬ 
netic fields tangential to the front and generate turbulence inside 
the front. This turbulence can re-accelerate the electrons and will 
result in a steeper radio index ( [ZuHone et al.|201^ . It’s not clear 
whether the sub-cluster hosted a radio-halo pre-merger. However, 
BCG 2, which is associated with the sub-cluster, is a radio galaxy 
(bottom panel Fig. [^, and may have provided the population of 
electrons to be re-accelerated by the sloshing of its core as it passes 
through the ICM of the main cluster. 

The survivability of a dense core under these conditions is not 
obvious. While the sloshing core and magnetic draping, at both the 
leading edge and due to the sloshing, may help confine the dense 
core gas and prevent much stripping and mixing, turbulence is also 
generated within the front which would act to disrupt the field and 
promote Kelvin-Helmholtz instabilities. 


5.4 inverse compton 

Inverse Compton (IC) emission from up-scattered CMB photons 
should be generated by the same population of relativistic elec¬ 
trons producing the radio halo. This non-thermal X-ray emission 
can in principle be detected as a power-law component or soft- 
excess at X-ray wavelengths (e.g. |Million & Allen|2009| l. The radio 
halo emission peaks in the core and the inferred high core X-ray 
temperatures could be indicative of an IC component to the spec¬ 
trum. Following |Million & Allen| ( |2009| ), we test this possibility 
by both allowing Nh to vary spatially, and by fitting an additional 
power-law component in regions with at least 15,000 counts. We 
find only one region, near the peak of X-ray emission, where Nh 
is significantly different from the Galactic value. If Nh is fixed to 
the Galactic value and an additional power-law component is fit, 
the best fit power law index requires an un-physical spectral index 
of a = —0.1 (the radio halo synchrotron index measured by |Orru| 
|et al.|2007| is a = 0.8 — 1.0). Whilst a two-temperature gas, with 
fixed abundance, does not offer a significant improvement in the 
goodness of fit to our data it is likely that the variation in Nh is 
caused by a multiphase gas associated with the dense core of the 
dominant subcluster. 


6 DISCUSSION 

The Abell 2219 system is undergoing a merger of two galaxy clus¬ 
ters in the NW-SE direction. The morphology of the X-ray emis¬ 
sion is elongated in this direction and three density discontinuities 
are observed along the projected merger axis. The system has a 
complicated core structure with multiple surface brightness edges 
and a long tail of stripped cooler gas to the SE. 

The most striking feature of the merger is a horseshoe of high 
temperature gas at 166 arcsec (595 kpc) towards the NW (Eig.[^, 
the leading edge of which coincides with a shock front (Eig.[7]and 
1^. In the reverse direction, at 155 arcsec (556 kpc), from the cen¬ 
troid, we detect, with >99% confidence in a temperature jump, a 
second shock front. The NW shock is associated with diffuse ra¬ 
dio emission, having a steeper spectral index than the remainder of 
the halo, yet the SE shock shows no evidence for associated radio 
emission. 

A probable cold front is observed in the SE wake of the 
stripped sub-cluster core. However, neither the surface brightness 
profile nor the deprojected density profile can reject at high con¬ 
fidence the null hypothesis that no jump exists (although the co¬ 


incidence of the feature with a radio edge adds evidence for the 
discontinuity being real, see Se ction | 5 .3|and Eig.[^. 

Recently, [Roediger et al.| ( |2014| , simulated the disruption of 
the atmospheres of elliptical galaxies falling into clusters. The au¬ 
thors found much of the gas, forming the wake, remains unmixed 
as a ‘remnant tail’ for long periods and, within this wake, sloshing 
cold fronts are expected due to the ‘push’ of the ICM. Behind the 
‘remnant tail’ is a deadwater region and stripping and mixing of 
gas only occurs beyond this deadwater region. Whilst these simu¬ 
lations were for individual galaxies, we may be observing a similar 
process in the stripped atmosphere of the sub-cluster. 

The long, low entropy and temperature wake in the SE sug¬ 
gests there has been little mixing of the tail gas with the surround¬ 
ing hotter ICM. The pressure map is smooth beyond the tail edges 
indicating pressure equilibrium is maintained with the surrounding 
gas. The feature observed at 359 kpc might be a sloshing cold front, 
within the remnant tail, or could indicate the end of the remnant tail, 
where the deadwater region occurs. 


6.1 Hot core and multiple surface brightness edges 


[Million et al.| ( |2010| ) identified a hot core in the A2219 system. Our 
observations confirm a temperature increase in several regions in 
the core. Notably, both BCGs, appear to be associated with high 
temperature regions, with the dominant BCG 1 associated with the 
highest temperatures (~15 keV), and densest gas (0.01 cm“^ ), in¬ 
dicating that the merger shocks are able to penetrate the central 
regions. 

A series of ‘inner’ surface brightness edges are also observed 
close to the centroid of the system (within 50” or ~ 180 kpc). Some 
edges are likely cold fronts from the sloshing cores (e.g. ‘bridges’ 
and ‘edges’ observed in the simulations of [Poole et al.|2006'] ) but a 
few of these edges are similarly associated with high temperature, 
high pressure regions and are likely driven by shocks around the 
core. We are unable to precisely constrain the deprojected temper¬ 
atures and densities in these regions due to the combination of the 
small widths of the shock regions and the high temperature gas. 
However, from the surface brightness discontinuities, we estimate 
Mach numbers of ~ 1.15. 

BCG 1 is not a strong radio source (see Eig.[9] bottom panel), 
and there are no indications in our X-ray surface brightness of 
‘holes’ resulting from bubbles blown by AGN jets, so the shocks 
are more likely a result of the dynamics of the merger rather than 
AGN driven shocks. The projected pressure maps show the dom¬ 
inant disruption occurring at the sides of the core, perpendicu¬ 
lar to the projected merger axis. This is where one might expect 
the shear fiow is largest and the gas less protected and more eas¬ 
ily stripp ed ([Roediger et al.[[2014[ ). However, from optical spec¬ 
troscopy, [Boschi^^r l[2004]r^ropose that the main A2219 core 
has also suffered the repeated in-fall of many small groups. In this 
scenario the surface brightness edges may be a consequence of 
many shocks experienced by the dense core due to the in-fall. 


6.2 Anatomy of a merger 

The X-ray morphology, the tail of cool gas and the outer surface 
brightness edges suggest the system is undergoing a significant 
merger event, with the main component along the NW-SE axis. 
This is consistent with the direction of elongation of the galaxy 
population. However, [Boschin et al.[ ( [2004[ ) also present evidence 
for a significant line-of-sight velocity gradient (~1000 km ) in 
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the E-W direction; from this they suggest the merger must have a 
component perpendicular to the plane of the sky, and that a ~45° 
angle is the most likely. The two BCGs also have a significant line- 
of-sight velocity difference (~ 1300 kms“^). If the merger is not 
exactly in the plane of the sky, projection effects would result in us 
under-estimating the shock Mach numbers. A higher Mach number 
would alleviate tensions between the radio edge emission and low 
Mach number NW shock. 

Simulations show that a forward propagating shock can have a 
‘horseshoe’ structure around core passage, with these features be¬ 
ing wiped out after ~1 Gyr ( [Ricker & Sarazin|2001[ [Poole et al.| 
[2006[[SEllman et al.|2008[ ). The strongest shock is observed lead¬ 
ing the less massive sub-cluster. Our projected temperature map 
shows a similar configuration, supporting the hypothesis that this 
is a very recent merger (~at core passage) - by non-equal mass - 
clusters, with the smaller cluster approaching from the SE (strong 
lensing also suggests the smaller sub-cluster is in the SE; [Smith[ 
[et al.[200^ . As it seems unlikely that the smaller subcluster would 
harbour the larger BCG the merger is probably being observed just 
before closest approach. At a projected distance of 595 kpc, the 
shock velocity translates to a time since core passage of the shock 
of ~ 0.26 Gyrs, supporting the scenario of a youthful merger. 


7 CONCLUSIONS 

Abell 2219 is a system in the throes of a violent merger, with the 
dominant component of the motion in the NW-SE direction. Both 
shock fronts and cold fronts in this system appear to be influencing 
the radio halo morphology and structure. The key results from our 
Chandra analysis are: 

(i) The projected temperature maps show that the cores are sur¬ 
rounded by a high temperature ‘horseshoe’. In the NW, this is co¬ 
incident with the surface brightness edge and density discontinuity 
which define the outer shock front. The shock front is coincident 
with a radio halo edge with a steeper spectrum than the remainder 
of the halo. The shock Mach number is M ~ 1.2, leading to a ve¬ 
locity of ~2100 km and time since core passage of the front of 
0.26 Gyr. 

(ii) A reverse shock with similar low Mach number is observed 
in the SE. This shock is not associated with any radio emission 
indicating that shock induced turbulence is not yet widespread in 
the cluster. 

(iii) Simulations of cluster mergers at core passage show similar 
‘horseshoe’ temperature structures and indicate that this is a early 
stage merger between non-equal mass ratio clusters. The youthful 
nature of the merger is supported by the positions of the BCGs; the 
first and second BCG are ~180 kpc apart and have probably yet to 
cross paths. 

(iv) We present the first evidence for a sloshing cold front ob¬ 
served in the ‘remnant tail’ of an in-falling sub-cluster. The cold 
front is observed in the SE and bounds the diffuse radio halo. The 
spectral index may steepen at this boundary indicating the cold 
front is amplifying the field at this location. 

(v) Both the X-ray cores associated with the BCG 1 and 2 are 
‘hot’ and multiple surface brightness edges surround them. These 
edges are likely of both shock and cold front origin and correspond 
to Mach numbers of ~ 1.15. 

The shock fronts in A2219 are not at great distances (the out¬ 
ermost shock being <600 kpc from the centroid and probable in¬ 
ner shocks within 180 kpc), so they have yet to propagate to the 


very low surface brightness regions of the system. As such A2219 
presents an interesting opportunity to study the detailed connec¬ 
tion of the thermal X-ray population and the non-thermal radio syn¬ 
chrotron emission. A population of relativistic electrons must exist 
in the system, generated by the three known radio galaxies, and the 
merging ICM is likely to be highly turbulent. However, the corre¬ 
spondence of the radio emission with the NW X-ray shock, and 
the spectral index variations in the core, where we see evidence 
for multiple shocks is tantalizing; low Mach number shocks may 
be more efficient at promoting synchrotron emission than currently 
thought. The tail of low temperature X-ray gas, sharp radio edge at 
the location of the cold front and lack of radio emission associated 
with the reverse shock suggest there is little mixing and turbulence 
behind the in-falling core. High resolution X-ray data combined 
with radio imaging and polarisation data should enable us to deter¬ 
mine the role the shocks and cold fronts are playing in compressing 
the cluster fields and re/accelerating the electrons and ions. 
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